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1 INTRODUCTION 



Since the concept of stellar populations was invented (|Baadel 
1 19441 '). stellar population synthesis (SPS) has been proven 
to be a useful tool for revealing clues about galaxy for- 
mation (|Tinslevlll96cill978l ; lBruzual fc CharlotJ 11993 ). Re- 
search on several crucial factors of stellar modelling, like con- 
vection, o pacity, heavy-elemen t mixing, helium content, and 
mass loss ( Chariot et al. 1996), resulted in increasingly accu- 
rate SPS models jBruzual fc CharlotJll993l; IWorthevlll994l; 
Bertelli et "all Il994l; iBruzual fc Chariot] 120031 ; iDotter et"al] 



20081 ; iBertelli et al l I2008I . 120091 ). In this scheme, estima- 
tion of single burst equivalent age and metallicity, was 
accomplished by comparing integrated light of observed 
galaxy with that of SPS models. However, most of the 
colours and absorption feature strength appear identical 
if the changes of age and metallicity satisfy S log (age) « 
—3/2 Slog(Z). This so-called age- metallicity degeneracy 
(Worthcy 1994) blocks the way to estimating galaxy ac- 
curate ages from optical broad band colours. At spec- 



ABSTRACT 

In the study of galaxy integrated light, if photometric indicators could extract age 
and metallicity information of high enough quality, photometry might be vastly 
more efficient than spectroscopy for the same astrophysical goals. Toward this end, 
we search three photometric systems: David Dunlap Observatory (DDO), Beijing- 
Arizona- Taiwan-Connecticut (BATC), and Stromgren systems for their ability to dis- 
entangle age and abundance effects. Only the Stromgren [ci] vs. [mi] plot shows 
moderate age- metallicity disentanglement. We also add to the discussion of optical 
to near-infrared Johnson- Cousins broad band colours, finding a great decrease in age 
sensitivity when updated isochrones are used. 
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in the optical so one spectrograph can cover them all) and 
the model grids in the observed index-index space open up 
to nearly orthogonal grids rather than collapsed into linear, 
degenerate overlapping segments. 

Spectroscopy is effective, but photometry also has its 
strong points. It is much more efficient for faint objects and 
low surface brightness galaxies. In the epoch of large sky 
surveys, like Sloan Digital Sky Survey (SDSS), Two Micron 
All Sky Survey (2MASS), Spitzer Space Telescope (SST), 
Large Synoptic Survey Telescope (LSST), and James Webb 
Space Telescope (JWST), photometry is undoubtedly more 
feasible than spectrometry. In fact, optical near-infrared 
(near-IR) colours have been exploited to possibly con- 
strain age and metallicity by several studies (|Peletier et al.l 
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or metallicity sensitive Lick indices ( 


— — (->- ~~ 
Burstein et al. 


1984; 


Faber et al.l 19851: Gorgas et al. 19931; 


Worthev et al. 


1994: 


Worthev & Ottavianil Il997l; iThomas et al.ll2003h. taken in 



pairs, are effective at breaking this degeneracy. Balmer in- 
dices vs. Fe-peak indices are widely used, for example, be- 
cause they are relatively easy to observe (all features are 
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1 199d;lde Jondl 19961 ; [Bell fc de Jone)l2000l ; ICarter et alj|2009[ 
Iconrov et al.ll2009l ). lc"ardiel et all (120031 ) suggested (V - K) 
colour shows g ood balance betwee n parameter degeneracy 
and sensitivity. IPessev et alj (2008) found that four age in- 
tervals (>10 Gyr, 2-9 Gyr, 1-2 Gyr, and 0.2-1 Gyr) are 
clearly separated in the optical near-IR colour-colour plot. 
These results hinge upon the models, and specifically upon 
the number of thermally pulsing asymptotic giant branch 
(TP-AGB) stars, which dominate near-IR light be tween 0.3 
and 2 Gyr of age (|Marastonll2"u05l ; iLee et al.li2007l ). 

The present work looks for age-sensitive or metallicity- 
sensitive indices in three ph otometric systems: David 
Dunlap Observatory (DDO, iMcClure fc Van den Berghl 
Il968l ; IMcClure! 1 19761) . Beijing- Arizona-Tai wan-Connecticut 



BATC, Fan et 



Stromgren 19661 : 



1996; Shan. 



Rakos et 



,ng 



et al.ll 19 98;). and Stromgren 
199ll .[l99d) systems. These 
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Table 1. Metallicity sensitivities at 8 Gyr, solar metallicity, for 
our best candidate pairs of colors 



Index 




DDO 0(41-42) 


-0.60 


DDO 0(35-48) 


1.42 


BATC 6075-6660 


1.15 


BATC 5795-6075 


1.97 


Stromgren [mi] 


1.23 


Stromgren [c±] 


2.41 



filters are narrower in wavelength than Johnson-Cousins 
broad band filters, but resemble the N IRC am narrow fil- 
terfl on board JWST. Can we find indices which are prac- 
tical to break the age-metallicity degeneracy? The tables 
and graphs in § 12. II explore this issue. In addition, equipped 
with the latest stellar evolut i onary isochrones of the Padova 
groujH (|Bertelli et all I20081 . I2009I . hereafter B09), we plot 
the model grids in (B — V) vs. (V — K) space, and com- 
pare with previous models in § 12.31 Observables from three 
samples are chosen to verify the credibility of parameters 
given by the model grids. Discrepancies between two model 
grids are briefly discussed in §[3] and a brief summary of the 
results and outlook for the future is given in §[3] 

2 ANALYSIS 

2.1 DDO, BATC & Stromgren 

In order to find the age-sensitive or metallicity-sensitive in- 
dex, we calculate the the metalli city-to-ag e sensitivity pa- 
rameter (Z ap ) using the evolving Worthevl (Il994l. hereafter 
W94) models, following IWorthevI (|l994l ) and IServen et all 
(|201ll ). The spectral library is synthetic and is sensitive 
to individual elemental abundances in the 300 to 1000 nm 
wavelength range. The metallicity sensitivity parameter is: 

_ 5I m /8\og(Z) 
sp 8I a / 5 log (age) 1 ' 

Where SI m /8\og (Z) is the metallicity partial derivative of 
the index at 8 Gyr, and 8I a /5 log (age) is the age partial 
derivative of the index at solar metallicity. Thus, Z 3p is the 
metallicity sensitivity at 8 Gyr, solar metallicity. 

Z 3p is surveyed for indices in three photometric sys- 
tem^. We refer the readers to related papers for de tailed de- 
script ion of filters and indices (See JTJ. As shown in lWorthevi 
|l994h . indices with Z ap close to 1.5 are less likely to break 
the age-metallicity degeneracy. Hence, indices with extreme 
Z sp values are selected as candidates and listed in Table 
What is not quantified in Table [1] is the dynamic range 
(compared to observational error) of a given index, and the 
dynamic range is low in DDO C(41-42), alas. 

1 http:/ /ircamera. as. arizona.edu/nircam/ 

2 http: //stev.oapd.inaf.it/YZVAR/ 

3 http:/ /astro. wsu.edu/models/isochrones. html 

4 Since the K band is suggested to be sensitive to young stel- 
lar populations, nonstandard colours relative to the K band, 
like DD041-K, BATC6075-K, were also examined. However, 
no promising index is found. 




C(35-48) 

Figure 1. DDO 0(35-48) vs. 0(41-42) plot. SSPs of the same age 
are connected as isochrones. Our models (W94) are given at age= 
1.5, 2, 3, 5 with log (Z) = -0.225, 0, 0.25, 0.50, and age= 8, 12, 17 
Gyr with log(Z) = -2.0,-1.5,-1.0,-0.50,-0.25,0,0.25,0.50. 
Solar metallicity SSPs are marked as filled squares to guide the 
eye. A vector for Ay = 1.0 mag is sketched to estimate the ex- 
tinction effect. 



• DDO system 

The extremely low Z sp value of C(41-42) colour makes it a 
promising age indicator to break the degeneracy. Figure [T] 
shows C(35-38) vs. C(41-42) plot of our models using W94 
stellar evolution. The isochrones are well separated, except 
SSPs with very low metallicities. However, the major prob- 
lem is small dynamic range («0.05 mag) of C(41-42). A 
little observational error or dust content (Ay- = 1.0 mag 
is sketched in Figure Q] to estimate the extinction effect) 
will destroy the well spaced isochrones. The short distance 
between the central wavelengths of DDO 41 and 42 filters 
(«100 A) is the reason for this defect. 

• BATC system 

BATC 6075-6660 and 5795-6075 colours show potential of 
breaking the degeneracy for their Z ap values. However, Fig- 
ure[2]does not work out as expected. Though the super-solar 
metallicity SSP isochrones are separated as indicated by Z ap , 
the sub-solar metallicity SSP isochrones are still highly de- 
generate. We note that the 6660 filter contains Ha, which 
explains age sensitivity of the 6075-6660 colour. However, 
we also note that the dynamic range of both colours is small 
compared to the Ay = 1.0 mag extinction vector, and that 
the 6075-6660 colour only spans 0.06 mag total, which seems 
small compared to reasonable observational errors. Further- 
more, since it is Hcv that is driving the age sensitivity of this 
colour, it will be sensitive also to nebular emission, which 
can easily overwhelm the stellar absorption in many galax- 
ies. 
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0.12 0.14 0.16 0.18 

5795-6075 



Figure 2. BATC 5795-6075 vs. 6075-6660 plot. Isochrones are 
described in Figure[T] A vector for = 1.0 mag is also sketched 
to estimate the extinction effect. 
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Figure 4. Stromgren c\ vs. mi plot of B09 models. The mod- 
els are given at age= 0.1,1,2,5,10,17 Gyr with log (Z) = 
-2.23, -1.23, -0.53, 0, 0.37. Solar metallicity SSPs are marked as 
filled squares to guide the eye. Observable symbols are the same 
as in Figure [3] 
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Figure 3. Stromgren ci vs. mi plot. Isochrones are described in 
Figure [l] Ages are not labelled for narrow space, but they can 
be easily inferred from Figure [4] Observables from three samples 
are l abelled with different sy mbols: Filled Triangles: Es and SOs 
fromlSc homber t et al] l|l993T ). Crosses: Nea rby Es and GCs fr om 
iRakos et alj Jl990T) . Filled Circles: Es from lOdell et al.l d2002h . 



• Stromgren system 

Since [ci] and [mijfjare reddening independent (jStromgren I 
Il966l ). and they also have pr omising Z sv values, we take a 
close look at these two indices. [Rakos et al.1 (|l99lh suggested 



ci 3 (u — v) — (v — b), mi = (v — b) — (b — y); 
[a] =ci- 0.2(6 - y), [mi] = mi + 0.18(6 - y). 



mi is a metal-line idex, while ci is a Balmer discontinuity in- 
dex. Figure |3] confirms [mi] as a metallicity sensitive index, 
because the SSP isochrones lie almost parallel to the [mi] 
axis. With its reasonable dynamic range (~0.6 mag) and 
its reddening free advantage, [ci] vs. [mi] plot is a promis- 
ing candidate index that might break the age-metallicity 
degeneracy. Three samples are selected to test the compe- 
tency of_ti*_JcjJ_j[S ; JmiJ_pjlot: (1) Ellipticals (Es) and SOs 
from ISchombert et all (|l993h . shown as filled Triangles in 
Figure [3] (2 ) Near by Es and globular clusters (GCs) from 
IRakos et al.l (|l990h . The y are labelled sep arately for clar- 
ity (crosses); (3) Es from lOdell et al.l (|2002l ). shown as filled 
circles. Note that the last two samples are recorded in a 
modified St romgren filter syst em called uz, vz, bz, yz. The 
relations of Ra kos et al.l (| 1996T) are employed to convert the 
observables to the uvby system. Figure [3] shows that our 
model grids based on W94 evolution do not fit the observa- 
tions well. B09 model grids are also plotted along the ob- 
servables in Figure [4] Comparing with Figure |3j the coverage 
improves greatly in the later model grids, but isochrones of 
age between 5 and 17 Gyr are still partially degenerate. Ob- 
viously, the [mi] indices of old metal-rich SSPs in B09 mod- 
els are systematically larger than that in W94 models. We 
also note that two sets of tracks are not the same, e.g., young 
age, low metallicity SSPs are not included in W94 models. 
To investigate the topic furthermore , we draw the [ci] vs. 
[mi] plot of B94 models (§ 12. 3[) . which looks strikingly sim- 
ilar to Figure [4] (not shown in this paper). Thus, the stellar 
phases causing the dramatic discrepancies between Figure [3] 
and|4]must be different in W94 and B09 models, but similar 
in B94 and B09 models. We defer discussion of the specific 
stellar phases because it is beyond the scope of this work. 
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Table 2. Index changes at 8 Gyr, solar metallicity 



Index 


C 


N 


O 


Na 


Mg 


Fc 


DDO C(41-42) 


0.028 


0.039 


-0.018 


-0.001 


-0.002 


-0.001 


DDO C(35-48) 


0.068 


0.026 


-0.033 


-0.005 


-0.002 


0.045 


BATC 6075-6660 


-0.009 


-0.006 


0.002 


0.000 


0.001 


0.002 


BATC 5795-6075 


0.005 


0.001 


-0.004 


0.004 


-0.001 


0.001 


Stromgren [mi] 


0.024 


0.021 


-0.025 


-0.003 


-0.032 


0.029 


Stromgren [ci] 


-0.078 


-0.015 


0.054 


0.003 


0.080 


-0.009 



2.2 Element Sensitivity 

In a non-solar scaled environment, indices may reflect those 
abundance changes. To trace the effects of different elements, 
we increase the stellar atmosphere model abundances of six 
major elements (C, N, O, Na, Mg, Fe) by 0.3 de^E one single 
element at a time. The sensitivity is modeled only in the 
stellar library; the isochrones are kept as published. Tabled 
outlines the index changes at 8 Gyr, solar metallicity. 

• The prominent carbon and nitrogen effect on DDO 
C(41-42) colour is not surprising, since the DDO 41 band 
overlaps CNi and CN2 Lick indices in wavelength. Fe and 
CN absorption lines are the main features under 4000 A, 
which explains the DDO C (35-48) colour changes at super- 
solar C, N, and Fe abundances. The behavior of increased 
O is opposite due to molecular balancing: The CO molecule 
has the strongest binding energy and so adding O decreases 
the C supply, weakening C2 and CN features. Thus, adding 
more O causes bluer DDO C(41-42) and C(35-38) colours. 

• The BATC 6075-6660 and 5795-6075 colours are rela- 
tively unaffected by six major elements, which makes them 
good colours to study systems with unknown element en- 
hancements. What we find matches the basic philosophy of 
creating this photometric system: avoiding known bright fea- 
ture lines and sky lines, in ord er to study the continuum of 
the spectrum (|Fan et al.lll996h . 

• Stromgren [mi] and [ci] are different from the four 
colours above, because they are the subtractions of two 
colours, [mi] is found to be metallicity sensitive in § 12.11 
which is consistent with its behavior in Table Five ele- 
ments have substantial influences on [mi]. The trend of the 
[ci] index change due to C and O abundance variations is 
reverse compared to other indices. Since u and v bands are 
both blanketed by CN absorption lines, our models predict 
that the impact of CN on the v band is greater than that 
on the u band. Therefore, higher CN abundance means a 
larger magnitude increase in the v band than that in the u 
band, leading to a smaller [ci] value. The strong Mg effect 
on [d] is echoed by th e Mg3835 (u band) and Mg4780 (b 
band) indices defined in lServen et all |2005l ). 



(l2007h inspected the ( B - V) vs. (V - K) plot us- 
ing iBruzual fe Chariot] (120031. hereafte r BC03), iMarastonl 
<|2005h . and BaSTI l|Cordier et all 120071 ) models. They sug- 
gested different treatments of convective core overshooting 
and differences in the TP-AGB phase cau se the inconsis- 
tency among different model grids (also see lMaraston et"ai1 
2006). To further discuss the robustness of the models, we in- 
vestigate the isochrones that looked so promising (B94) and 
a later, even more complete version (B09) from the same 
group. 

• B94 models 

Perhaps partly because of the inclusion of all stellar evolu- 
tionary phases, B94 models were widely adopted as basic 
stellar evolution isochrones for e volutionary synthesis mod- 
els, s uch as: BC03; PEGASE dFioc fe Rocca-Volmerange] 
1997 1 . Il999l ); and STARBURST99 i|Vazquez fe Leithererl 
20051) models. Figure [5] shows the (B - V) vs. (V - 



K) plot of B94 models. Three observational samples are 
selected for c o mpar ison: (1) old, metal-poor GCs from 
iBurstein et alj (1 19841) (filled circles); (2) old, metal-rich Es 
from lPeletierl l|l989l ) (filled triangles); (3) SWrfl type IV- VI 
Large Magellanic Cloud (LMC) star clusters (crosses). The 
reddening-corre c ted (B — V) and (V — K) colours come from 
IVan den Berghl l|l98ll) and IPersson et all l| 19831 ) . The last 
sample is selected a s "intermediat e age" SSPs, roughly 1 to 
8 Gyr according to IPersson et alj (|l983l ). Generally speak- 
ing, observables of these three samples locate at expected 
locations on the model grids. 
• B09 models 

With improved understanding of TP-AGB phase, Padova 
group updated the stellar evolutio n models. Compared with 
previous models of this group (B94: lGirardi et al.l2000l ). B09 
models adopted a sophisticated TP-AGB model which in- 
cludes third dredge -up, hot bottom bu rning, and variable 
molecular opacities (|Marigo et alj |2008). Also, basic observ- 
ables of AGB stars in the Ma gellanic Clouds (MCs) we re re- 
produced by the new models l|Girardi fc Marigd[2007l ). The 
(B - V) vs. (V - K) plot of the B09 models are shown 
in Figure [6] Observables are also plotted as in Figure [5] 
Isochrones of age larger than 5 Gyr are consistent with that 
of the B94 models. Es and GCs locate on the metal-rich 
and metal-poor sides, respectively. However, isochrones of 
1 and 2 Gyr change substantially: the (V — K) colour of 
log (Z) — —0.53 becomes redder than that of solar metal- 
licity. An unreported "S" shape is found if log (Z) — —0.53 
SSPs are connected. Another unexpected change is about 
the isochrone of the youngest age, 100 Myr: the (V — K) 
colour of the lowest metallicity (log (Z) = —2.23) is even 
redder that of the highest metallicity (log (Z) = 0.37). 



2.3 (B - V") vs. (V - K) Plots 

The possibility of breaking age - metallic i ty de g eneracy with 
optical plus near-IR co lours jde Jond 1 19961 : | James et al.l 
120061 : iPessev et al.l 12008). and the upcoming era of near- 
IR telescopes (eg: JWST) inspire us to delve t he optical 
near-IR colour— colour plot. Prior to this work, iLee et all 

6 0.15 dex for C, to avoid the danger of turning to a carbon star 
llServen et al.ll2005l) . 



3 DISCUSSION 

3.1 Stromgren System: Empirical Leverage? 

Prior to this work, c\ or mi were paired wi th (u — v), (b — y), 
or other colours to estim ate galaxy age (|Russell fe Bessell 
Il989l: iRakos et alj 1 1996') . In this paper, the extinction in- 
dependent Q-indices [ci] and [mi] are assembled, but do a 

7 ISearle et al.l dl980l) 
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Figure 5. (B - V) vs. (V - K) plot of B94 models. The 
models are given at age= 0.1,1,2,5,10,17 Gyr with log (Z) = 
— 1.7, —0.7, —0.4, 0, 0.4. Solar metallicity SSPs are marked as filled 
squares to guide the eye. SSPs of the same age are connected 
by solid lines, while SSPs of the same metallicity are connected 
by dotted lines. A vector for Ay = 1.0 mag is sketched to es- 
timate the extinction effect. Observables from three sa mples are 
labelle d by different symbols: Fi lled Triangles: Es fro m iPeletierl 
jl98Sh . Filled Circles: GCs from lBurstein et alj il984Tl . Crosses: 
SWB type IV— VI L arge Magellanic Cloud s tar clusters from 
I Van den Bergbl Jl98lf> and lPersson et alj Jl983f) . 
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Figure 6. [B - V) vs. (V - K) plot of B09 models. The 
models are given at age= 0.1,1,2,5,10,17 Gyr with log (Z) = 
—2.23,-1.23,-0.53,0,0.37. Other symbols are the same as in 
Figure [5] 
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Figure 7. Stromgren [ci] vs. [mi] plot of B09 models. The mod- 
els are discribed in Figure [4] Solar metallicity SSPs are marked 
as filled squares to guide the eye. Filled Circles: S and S+ of 
lOdell et alj d2002ft . 



mediocre job breaking the age-metallicity degeneracy. As 
seen in Figure |4j several observables lie outside the B09 
model grids. Three reasons are hypothesized for these out- 
liers: First, the conversion rel ations between uz, y z, bz, yz 
and uvby system presented in iRakos et all §996) are em- 
pirical, which inevitably increases the uncertainty; Second, 
the dynamic range of [ci] index is not large compared to the 
observational error; Third, isochrones of age between 5 and 
17 Gyr are partially degenerate, which reduces the span of 
model grids. 

The model isochrone of 100 Myr lies low i n the \ci] vs. 
[mi] p lot. This is similar to the statement in IRakos et al. 
( 1996) that most of the starburst galaxies have mz valued 
less than —0.2, and they are well separated from other types 
of galaxies. To test the potential of [mi] index as an star 
formation rate indicator, S ("star formation rates equivalent 
to normal disk ga laxy") and S+( "starburst objects") from 
lOdell et alj (j2002r i are selected and shown in Figure [7] Most 
of the mz values of this sample are less than —0.1, supporting 
the idea that starburst galaxies are bluer in [mi]. However, 
we detect no clear [mi] value that separates Es and star 
forming galaxies. 

Stromgren system photometry is far less developed than 
broad band photometry for galaxies, however. That, plus 
the large apparent observational scatter in Figures [3J U and 
[7]read us to conjecture that Stromgren photometry may still 
have great potential for breaking the age-metallicity degen- 
eracy, but higher quality data are needed to be certain. 



8 mz = (vz — bz) — (bz — yz) 
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3.2 Model Grid Differences 

Three SPS models are assembled in this work: W94, B94, 
and B09 (i.e., primitive, competent, and sophisticated, re- 
spectively.). Figure [3] [4] [5] and[6]clearly illustrate the model 
differences in [ci] vs. [mj and (B — V) vs. (V — K) plots. 
During the whole work, we calculate the magnitudes of 
each bands in three SPS models with the same codes, thus 
W94, B94, B09 models share the same spectral library, IMF 
<|Sarpeterlll955T ) with lower mass limit M mi „ = O.IMq, upper 
mass limit M max = 100M Q , and etc. Most of the model grid 
discrepancies rely on different stellar evolution descriptions: 

• W94 models are a merging of isochrones and tracks from 
many sources. Compared to B94, W94 models have a fac- 
tor of two more stars in the upper RGB. The core- helium 
burning phase was approximated by a single red clump. The 
luminosities and lifetimes of AGB stars were drawn from 
theoretical tracks. 

• B94 models included all phases of stellar evolution until 
the remnant stage in all mass ranges. Convective core over- 
shooting was also considered in the models. Nevertheless, 
the an alytic prescription of TP-AGB phase was approxi- 
mate. (|Marigo et alj|2008h 

• B09 models a dopted the latest TP-AGB model of 
iMarigo et alj j2008). Compared with previous ones, this TP- 
AGB model has included several crucial effects and it was 
calibrated by AGB stars in the MCs. Bertelli et al. reckoned 
B09 models as satisfactory achievement. 

In Figure [6] SSPs with log (Z) = —0.53 have the red- 
dest (V — K) colours along isochrones of 1 and 2 Gyr. This 
feature is not seen in B94 model grids. Since B94 and B09 
models mainly differ in TP-AGB model treatments, which 
affect the near- IR colour seriously between 0.3 and 2 Gyr 
l |Marastonl i2005). it is suggested different treatments of TP- 
AGB stars are responsible for the changes of model grids. 
Considering the metallicities of star clusters in the LMC 
(mean log [Z) = —0.42 for clusters with 1 ^ t < 2 Gyr 
and —0.57 for clusters with 2 ^ t < 10 Gyr, IPessev et al.l 
120081 ). SWB IV- VI star clusters are better fitted in B09 
models — B94 models give much higher metallicities. Nev- 
ertheless, noting that B09 models are calibrated by AGB 
stars from the MCs, and that the metallicities of star clus- 
ters in the LMC are around log (Z) = —0.53, it is suspected 
that SSPs of other metallicities, especially the extremely low 
metallicity, may not be well calibrated. The scarcity of very 
low metallicity AGB stars around solar neighbourhood im- 
pedes empirical calibration. It is reasonable to suggest closer 
scrutiny for SSPs with extremely low metallicities around 1 
and 2 Gyr. 

Another unexpected feature in Figure [S] is the twist of 
the 100 Myr isochrone. This youngest isochrone is unique: 
the (V — K) colour at the lowest metallicity is redder than 
that of the highest metallicity. This is inconsistent with the 
usual sense that metal-rich SSPs are redder than metal-poor 
SSPs (See the 100 Myr isochrone of Figure [5]). The change 
of the 100 Myr isochrones between B94 and B09 models 
challenge the method of estimating ages and metallicities 
of young SSPs (0.1 < t < 1 Gyr) from optical near-IR 
colour— colour plots (|Brocato et al. I ll999l : iHunt et alj|2003h . 

Finally, we check B94 and B09 model differences in 
the popular (B - R) vs. (R - K) plot jde Jongj Il99fj : 




B94 - 



0.5 1 1.5 2 

B-R 

Figure 8. [B - R) vs. (R - K) plot of B94 and B09 models. The 
B94 models are given at age= 0.1, 1, 2, 5, 10, 17 Gyr with log (Z) = 
— 1.7,-0.7,-0.4,0,0.4, while the B09 models are given at age= 
0.1, 1, 2, 5, 10, 17 Gyr with log (Z) = -2.23, -1.23, -0.53, 0, 0.37. 
Note that the two sets of tracks are not identical. Isochrones of 
B09 models are labelled by solid lines, and isochrones of B94 
models are labelled by dashed lines. Solar metallicity SSPs are 
marked as filled squares to guide the eye. A vector for = 1.0 
mag is sketched to estimate the extinction effect. 

iBell fc de Jo"^3 l200d . Figure |SJ). Solid lines are isochrones 
of B09 models; dashed lines are isochrones of B94 models. 
Differences between two model grids are similar to (B — V) 
vs. (V - K) plots: SSPs of log (Z) = -0.53 have the red- 
dest (R — K) colours along the 1 and 2 Gyr isochrones in 
B09 models; The (R — K) colour at the lowest metallicity 
is redder than that at the highest metallicity along the 100 
Myr isochrone of B09. Thus, these features are dependent 
upon near-IR band (K), instead of visible bands (V or R). 
Considering the TP-AGB treament differences of B94 and 
B09 models and the TP-AGB effects on K band, it is again 
suggested TP-AGB stars are responsible for the changes of 
youn g SSP isochrone s (0.1 < t < 2 Gyr), which is consistent 
with lMaraston et alj (|2006t ) and lLee et ail (|2007f ) 



4 SUMMARY 

Motivated by the need of estimating ages and metallicities 
from photometric systems, this work explored the DDO, 
BATC, Stromgren systems in the hope of finding age sen- 
sitive or metallicity sensitive indices. Three index— index 
plots were examined but only Stromgren [ci] vs. [mi] plot 
showed moderate age-metallicity separation. Four indices 
of the DDO and Stromgren systems changed substantially 
while increasing C and N abundances. O traces an oppo- 
site trend compared to C and N, since O tends to form 
CO with C. Then, we turned to the literature-posited age- 
metallicity disentangling space — the optical plus near-IR 
colour— colour plot. Updated stellar evolution gives no sup- 
port for clear disentanglement of age. Even mean age, there- 
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fore, is a subtle effect in colors, and would require, seem- 
ingly, both finer-tuned models and observations of greater 
accuracy. 
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